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ABSTRACT 

We stochastically formulate the theory of scalar quantum electrodynamics 
on a de Sitter background. This reproduces the leading infrared logarithms 
at each loop order. It also allows one to sum the series of leading infrared 
logarithms to obtain explicit, nonperturbative results about the late time 
behavior of the system. One consequence is confirmation of the conjecture 
by Davis, Dimopoulos, Prokopec and Tornkvist that super-horizon photons 
acquire mass during inflation. We compute ~ 3.2991 x H'^. The scalar 
stays perturbatively light with M"^ ~ 0.8961 x 3e^if^/87r^. Interestingly, the 
induced change in the cosmological constant is negative, (5A ~ —0.6551 x 
2>GH^/'K. 
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1 Introduction 



Gravitons and massless, minimally coupled (MMC) scalars are unique in be- 
ing massless without classical conformal invariance. The combination of these 
properties causes the accelerated expansion of spacetime during inflation to 
tear long wavelength virtual quanta out of the vacuum [U |2[ [3]. As more 
and more gravitons and MMC scalars emerge from the vacuum, the metric 
and MMC scalar field strengths experience a slow growth. The effect can be 
felt by any quantum field theory which involves either the undifferentiated 
metric or an undifferentiated MMC scalar. 

A typical example is afforded by the MMC scalar with a quartic self- 
interaction, 

£ = + 5Z)d,vd,vg^''^ - ^6^ip'Ry^ ~ l! + 5A)/v^ . (1) 

Consider this theory quantized on a nondynamical, locally de Sitter back- 
ground, 

ds"^ = -dt^ + a^dx ■ dx where a{t) = e^^ . (2) 

If the finite parts of the renormalization constants (and the cosmological 
counterterm) are chosen to make the expectation value of the stress tensor 
vanish at t = 0, then an explicit two loop computation using dimensional 
regularization reveals the following results for the induced energy density 
p{t) and pressure p{t) JH [5], 

Pit) = ^.iWia)]+0{X'), (3) 



(2vr) 



Pit) = (^{-^i^'H-i^M«)} + o(A2)- (4) 

The factors of ln(a) = Ht in expressions (jSMj) are known as infrared 
logarithms. They derive from the slow growth of the scalar field amplitude 
that is apparent even in the free theory P [3, [8], 

Qq (f'^ix) f2o) = Divergent Constant + - — -rln(a) . (5) 

The two loop expectation value of the stress tensor acquires two such factors 
coming from the —^)\^p^g^i, term. 
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Any quantum field theory which involves undifferentiated MMC scalars or 
metrics will show similar infrared logarithms in some of its Green's functions. 
They arise at one and two loop orders in the scalar self-mass-squared of this 
same theory [9],|10]. In scalar quantum electrodynamics they have been seen 
in the one loop vacuum polarization [HI [12] and the two loop expectation 
values of certain scalar bilinears [13] . In Yukawa theory they show up in the 
one loop fermion self-energy [T31 [TJ] and in the two loop coincident vertex 
function [16]. In pure quantum gravity they occur in the one loop graviton 
self-energy [T7] and in the two loop expectation value of the metric \18^. When 
quantum gravity is coupled to a massless, Dirac fermion they occur in the one 
loop fermion self-energy [TSUSD]- They even contaminate loop corrections to 
the power spectrum of cosmological perturbations [2T], [221 [221 [211 [2S1 [2S] and 
other fixed-momentum correlators |27j . 

Infrared logarithms introduce a fascinating secular element into the usual, 
static results of quantum field theory. For example, without infrared loga- 
rithms, the expectation value of the stress tensor of ([1]) would be a constant 
times g^i,. With the same renormalization conventions [H [S] in computing 
([3}[1D, the constant would actually be zero! 

The most intriguing property of infrared logarithms is their ability to 
compensate for powers of the loop counting parameter which suppress quan- 
tum loop effects. Indeed, the continued growth of ln(a) = Ht must eventually 
overwhelm the loop counting parameter, no matter how small it is. However, 
this does not necessarily mean that quantum loop effects become strong. 
The correct conclusion is rather that perturbation theory breaks down past 
a given point in time. One must employ a nonperturbative technique to 
follow what happens at later times. 

Certain models lend themselves to resummation schemes such as the 1/A^ 
expansion [2Sl 121] but a more general technique is suggested by the form of 
the expansion for p(t) in ([1]), 



oo 



£=2 



2t-2 



=i/^5]A^-^ Q,o ln(a) ln(a) + . . . + c,,2i-2\n\a)\. (6) 



2^-3 



Here the constants q ^ are pure numbers which are assumed to be of order 
one. The term in (jS]) involving [Aln^(a)]^~^ is the leading logarithm contri- 
bution at £ loop order; the other terms are subdominant logarithms. Pertur- 
bation theory breaks down when ln(a) ~ 1/V^, at which point the leading 
infrared logarithms at each loop order contribute numbers of order one times 
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H^. In contrast, the subleading logarithms are all suppressed by at least one 
factor of the small parameter ^1. So it makes sense to retain only the 
leading infrared logarithms, 



p(t) -^i/^ 5: Q,o [A ln^(a)j'"\ (7) 
e=2 

This is known as the leading logarithm approximation. 

Starobinskii has developed a simple stochastic formalism [30] which repro- 
duces the leading infrared logarithms at each order for any scalar potential 
model of the form, 

£ = _!(! + SZ)d^^d,^g>^''y^ - V{^),/^ . (8) 

Probabilistic representations of inflationary cosmology have been much stud- 
ied in order to understand initial conditions [31^^ r32| and global structure 
[331 131] ■ However, we wish here to focus on Starobinskii's technique as a 
wonderfully simple way of recovering the most important secular effects of 
inflationary quantum field theory [351 EHl [371 IM]- It is of particular impor- 
tance for us that Starobinskii and Yokoyama have shown how to take the 
late time limit of the series of leading infrared logarithms whenever the po- 
tential V{(p) is bounded below [32] • This is the true analogue of what the 
renormalization group accomplishes in fiat space quantum field theory and 
statistical mechanics. 

The solution of Starobinskii and Yokoyama is an amazing achievement, 
but it only gives us nonperturbative control over the infrared logarithms 
which arise in scalar potential models (jH]). The most general models which 
show infrared logarithms possess two complicating features: 

• Couplings to fields other than MMC scalars and gravitons; and 

• Interactions which involve differentiated MMC scalars and gravitonsQ 

An important step forward was a recent leading log solution for the model 
comprised by a MMC scalar which is Yukawa-coupled to a massless, Dirac 



"'^Of course there would be no infrared logarithms if all the MMC scalars and gravitons 
were differentiated. However, infrared logarithms must arise, in the expectation values of 
some operators, from interactions which involve at least one undifferentiated MMC scalar 
or graviton. Examples include the h^dhdh interaction of pure quantum gravity [171 118] 
and scalar interactions of the form Lp^dipdtp [2T1 141j. 
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fermion [TB] . That model possesses the first comphcating feature but not the 
second. In this paper we derive a similar leading log solution for MMC scalar 
quantum electrodynamics (SQED), 

£ = -{l + 5Z2){d^-ieA^)(t)*{d, + ieA,)(t)g^''^- 5i<P*<pR^ 

-^-{l + 6Z,)F,,F,,g^''g'^^^ - ^(0*0)2^ , (9) 

Although this model has derivative interactions we will see that a gauge 
choice permits one to avoid them at leading log order. We still do not have 
a full understanding of how to treat derivative interactions. 

Section 2 of this paper summarizes a recent all-order derivation [iQl HI] of 
the Starobinskii formalism for scalar potential models. Of course it is crucial 
to understand why the technique works in order to apply it to more general 
theories. That problem is discussed in section 3, reaching the conclusion that 
one deals with other fields by integrating them out and then stochastically 
simplifying the resulting effective action. Section 4 accomplishes this for the 
vector potential of SQED. One must also integrate out the vector potential 
from any operator whose expectation value is desired. The resulting, purely 
scalar operator is then stochastically simplified before computing the lead- 
ing log contribution to its expectation value using Starobinskii's technique. 
We do this in section 5 for the various constituents of the SQED stress ten- 
sor. Because the leading logarithm limit of SQED gives a model of the form 
(IH]), with a potential which is bounded below, one can exploit the solution 
of Starobinskii and Yokoyama to make explicit, nonperturbative predictions. 
We do this in section 6 for the expectation values of ip*{x)ip{x), F^y{x)Fp^{x) 
and T^y{x). We also confirm the remarkable conjecture of Davis, Dimopou- 
los, Prokopec and Tornkvist that super-horizon photons acquire mass during 
infiation |32l HS] . The eventual, nonperturbative photon mass-squared turns 
out to be about a hundred times larger than perturbative estimates. 



2 Deriving Starobinskii's Formalism 

Infrared logarithms arise in explicit perturbative computations of Green's 
functions formed from quantum field operators. It is an amazing fact that 
one can reproduce the leading infrared logarithms in any model of the form 
(IH]) using a formalism in which the fields are classical random variables. One 
wonders, what became of the Uncertainty Principle? What became of the 
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ultraviolet divergences and the counterterms? And how did the stochastic 
jitter emerge? 

In order to apply Starobinskii's formalism to more general models one 
must understand why it works. That is the task of this section. We begin by 
expressing the dimensionally regulated, Heisenberg field equations in Yang- 
Feldman form We then explain how infrared logarithms arise and, of 
crucial importance, the conditions for any expectation value of undifferen- 
tiated fields to receive a leading logarithm contribution. Based upon this 
understanding, we identify a series of simplifications that can be made to the 
Yang-Feldman equation without in any way affecting the leading infrared log- 
arithms. Taking the time derivative of this simplified Yang-Feldman equation 
results in Starobinskii's Langevin equation, with the white noise emerging as 
the time derivative of the simplified free field. The section closes with a re- 
view of the nonperturbative solution of Starobinskii and Yokoyama obtained 
[39] for the late time limit of any model with a potential which is bounded 
below. 

2.1 The Free Field Expansion 

The dimensionally regulated, Heisenberg operator equation for is, 

if + iD- l)Hcp -%ip + = . (10) 

Integrating it results in the Yang-Feldman equation [Hj, 

rt r V'(ip(x')) 

ipit,x) = Mt,x) - / dfa'""-' / rf^'VG(a:;xO . (11) 

JO J 1 + oZ 

A number of quantities in fill I) require definition. The free field ipo{x) and 
the retarded Green's function G{x; x') are, 

r d^~^h 

Mt,x) ^ j-^-^^e{k-H)[u{t,k)e^-'a{k) + u\t,k^^^^ (12) 

G{x- x') = ie{At)J-^^^j^e"'-^^{u{t, k)u*{t', k) - u*{t, k)u{t', k)} . (13) 
Here At = t — t', Ax = x — x' and the mode function u{t, k) is. 
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The nonzero commutation relations of the canonically normahzed creation 
and annihilation operators are, 



It follows that the retarded Green's function can be expressed as the com- 
mutator of two free fields n 



Iterating the Yang-Feldman equation generates the usual interaction pic- 
ture expansion of the field, in this case expressed in terms of a free field (foix) 
which agrees with the full field and its first time derivative at t = 0. Without 
worrying about operator ordering, we can write out the first few terms of this 
expansion. 



A nice diagrammatic representation for this expansion has recently been 
given by Musso [47] . 

The integrals over x'^ and x"^ are known as vertex integrations. The 
expectation value of any operator which involves <f{x) — for example, (p^{x) 
or ip{x)ip{x') — can obviously be reduced to a sum of terms, each one of 
which consists of a number of vertex integrations of Green's functions times 
the expectation value of some number of free fields. The expectation value of 
free fields can be further reduced to a sum of products of expectation values 
of two free fields. 

•^The restriction to fc = ||fc|| > _ff in the free field mode sum is imposed to avoid an 
infrared singularity in the free propagator [45j . The physical reason for this singularity is 
that no causal process would allow an experimenter to prepare the initial state in coherent 
Bunch-Davies vacuum over an infinite spatial section. Sensible physics can be regained 
either by employing an initial state for which the super-horizon modes are less strongly 
correlated [33], or else by working on a compact spatial manifold such as T^^^ for which 
there are initially no super- horizon modes [Tj. In both cases the modes with k < H are 
effectively absent. Note also that one typically removes the cutoff on any mode sum, such 
as ([T51) . which is not singular at A: = 0. 




(15) 



(16) 
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2.2 Genesis of Infrared Logarithms 



Infrared logarithms derive from two sources: the expectation values of pairs 
of free fields and the vertex integrations. 

• Free Field Expectation Values: Let us first consider the expectation 
value of two free fields in the presence of the state which obeys a{k)\VL) = 0, 



n 



n 



H 



B{k - H)e'' '^Mt, k)u*{t\ k) , (18) 

JD-■^{k^x)u{t,k)u*{t'k) 

>-• (19) 



dk k 



D~2 



)D-2 



TC 2 



At high k the Bessel functions in (1191) oscillate for ^ x'^, which makes 
the integral converge. Even at coincidence there is no possibility of an a- 
dependent ultraviolet divergence, 



(pl{x) 



)HaD-^ Jh 



dk k 



D-2 



^Ha^ 



H 



D-2 



2^7r— r 



dz z 



D-2 



H'iliz) 

2 



(20) 
(21) 



The only possible a dependence is the finite contribution from the infrared. 
The small k expansion of the integrand in (fT9l) is, 



k 



^__2 ikAx)u{t, k)u*{t', k) 



2D-2^^(k^)^ 



■D-l 



H 



D-2 



2n°^^ k 



O 



k' 



ma^ ' ma'^ ' 



k^/\x' 



(22) 



Had the lower limit not been cut o?i aX k = H this would give a logarithmic 
divergence. With the infrared cutoff there is no divergence, but one does get 
a large logarithm. It derives exclusively from the first term of fl22l) , integrated 
up to the point where the expansion breaks down and the integrand begins 
to oscillate. This point is /c ~ Min(ifa, ifa', 1/Ax). By taking account of 
causality. 



Ax < 



1 



1 



Ha Ha' 



> Mm{Ha, Ha') 



(23) 
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we see that the upper hmit is actually k = Min{Ha, Ha') 



r(^)if^-2 ^Hadk 



Itt 2 



T{D-1] 



H k (47r)f r(f) 



21n(a) 



47r^ 



ln(a) (in D = A) . 

(24) 

To summarize, infrared logarithms from ipQ{t,x) derive exclusively from the 
range H < k < Ha{t). Further, only the first term in the long wavelength 
expansion of the mode functions contributes. 

• Vertex Integrations: We turn now to the second source of infrared 
logarithms, which is vertex integrations. It is apparent from expressions (fT3l - 
H]) that the most infrared-singular part of the mode function drops out of 
the retarded Green's function, 



ie{At) [u{t, k)u*{t', k) - M*(t, k)u{t', k) 

e{At) r 1 



{D-1)H 



-.rD-l 



,D-1 



1+0 



k^ 



(25) 



Hence the retarded Green's function cannot contribute infrared logarithms. 
However, consider the vertex integration of G{x; x') against n powers of In(a'), 



J^dt'a""-^ J d''-^x'G{x;x')[\n{a') 



1 



{D - 1)H Jo 
1 

(n 



fdt' fln(a') 



(26) 



1)(D- l)iJ2 



ln(a) 



O [ln(a)]" . (27) 



The temporal vertex integration has increased the number of infrared loga- 
rithms from n to n + 1. 

A temporal vertex integration can only produce an additional infrared 
logarithm when it receives nearly equally weighted contributions from its full 
range. This requires that no factors of a' should remain after multiplying by 



■^We have also exploited the doubling formula [48 to write, 

'L>-1\ ^ r(D-i) 



r 



2^-2 r(f) 
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the a'^ ^ from the measure and performing the spatial vertex integrations. 
For example, consider the two terms from the right hand side of f l26l) . 

f dt' \\n{a')f 
Jo 

[' dt' (-)'''\ln{a')r 
Jo ^ 



To summarize, infrared logarithms from vertex integrations derive entirely 
from the part of the long wavelength expansion of the Green's function which 
goes like 1/a'^^^. 

2.3 Conditions for a Leading Log Contribution 

We have seen the ways in which infrared logarithms originate from the ex- 
pectation values of pairs of free fields, and from vertex integrations. It is 
now time to consider the crucial issue of how many infrared logarithms must 
derive from each source in order to reach leading log order. Although the 
answer is completely general, it is easier to explain in the context of the 
quartic self-interaction, 

V{^) = ^6^^'R+^{X + 6X)^'. (31) 

Let us first establish that the various counterterms cannot make leading 
log contributions. This follows from the number of fields they carry and the 
number of factors of A they involve in this model [H [5], [9], 

6Z = 0(A) , 6^ = 0(A) , 6X = 0{X^) . (32) 

The coupling constant renormalization, 6X, multiplies the same cp^ term as 
A, so contributions involving it produce the same structure of infrared log- 
arithms as contributions from Ayj^. Because contributions from 6X(p^ have 
at least one extra factor of A, with no more infrared logarithms, they can 
never be leading order. The same argument applies to the field strength 
renormalization, 6Z, on account of the fact that it enters the field equations 
in the form, V'{!f)/{1 + 6Z). While both 6C, and A go like A, the former 



Ma)] 



n+l 



dt' {Ht'Y 

[ln(a)]" 
(D - l)H ~ (D - IfH 



{n + l)H 

e^^* - 1 " 

n[ln(a)]"-i 



+ . . . . 



(28) 
(29) 

(30) 
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multiplies (p^ while the latter multiplies ip^. So no contribution involving 5^ 
can produce as many infrared logarithms, at the same order in A, as con- 
tributions involving only the Ay?^ term. We therefore get exactly the same 
leading log contributions from the simplified Yang-Feldman equation without 
the counterterms, 

^{x) = (^o(x) - - r dt' a""^' / d""-^^' G{x- x')^\x') . (33) 
6 JO J 

Let us consider the generic form of the free field expansion derived from 
( l33|) . keeping track only of the number of A's, the number of vertex integra- 
tions, and the number of fields at any point. We already have symbols for 
the coupling constant A and the field v?o- Let us employ the symbol "J" to 
denote generic vertex integrations, 

/ = ^* dt' a'^-^ j d^'-^x'Gix; x') . (34) 

In this notation we might render (l33l) as follows, 

99 ~ (yjQ + XI^P^ ■ (35) 

Note that we do not worry about signs or numerical factors such as 1/6, nor 
do we worry about which spacetime points the various fields reside. In this 
generic language it is simple to iterate fl35l) to exhibit the generic form of the 
free field expansion, 

~ (^0 + XIvl + A'/Vo + A^^Vo + • • • • (36) 

In other words, each additional factor of A involves one more vertex integra- 
tion and two more free fields. 

The same generic form (136!) applies to any operator whose VEV we might 
wish to compute. For example, the product of N <^'s — even at all different 
points — would be rendered, 

~ (v.o)^{l + M^l + ^'iVo + X'lVo + •■•}• (37) 

Note that we do not worry, at this level, about which spacetime points the 
various free fields reside, or which of the parent full fields gives rise to correc- 
tions in the free field expansion. The only interesting things are the number 
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of factors of A, the number of vertex integrations, and the number of free 
fields. 

Now recall that leading log corrections to any result in this theory must 
produce two infrared logarithms for each extra factor of the loop-counting 
parameter A [H [5], [9l [10]. From the generic free field expansion fl371) we see 
that each additional factor of A is accompanied by one new vertex integra- 
tion and two new free fields. Because a vertex integration can add at most a 
single infrared logarithm, as can any pair of free fields, it follows that leading 
log contributions to the expectation value of any combination of undifferen- 
tiated full fields (fi can only result when each vertex integration and each free 
field contributes to an infrared logarithm. That turns out to be a general 
conclusion for any potential V{(f). 



2.4 Starobinskii's Langevin Equation 

It is well to summarize what we have learned in the previous sub-sections: 

1. For the expectation value of any combination of undifferentiated full 
fields to receive a leading log contribution, every {ipo{x)(po{x')) and 
every vertex integration must contribute an infrared logarithm. 

2. The infrared logarithm from ipo{t, x) derives exclusively from modes in 
the range H < k < Ha{t), and from just the leading term in the small 
k expansion of u{t,k). That is, the following replacement makes no 
change at leading log order, 

d^'-'k 



ip.it, — / ^^-^ e{k - H)e{Ha - k) 



X 



^A-nH ^ k 



— ) ' {e'^-"^a(fc) + e-'^-V(A?)}. (38) 



3. The infrared logarithm from a vertex integration derives exclusively 
from the term in the small k expansion of u{t, k)u*{t', k)—u*{t, k)u{t', k) 
which goes like l/a'^~^. That is, the following replacement makes no 
change at leading log order. 
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Because of point 1 we can make the simplifications fl38ti39p in the Yang- 
Feldman equation. Because the infrared truncation of (1381) removes any pos- 
sibihty for uhraviolet divergences, we can set D = 4 to write the infrared 
truncated free field as, 

r d^k H - 

$o(t,x)= y ^^(A;-i/Mi/a-A:)-^{e^'=-a(^) + e-^-V(^)}. (40) 

Note that we have used a new symbol, ^q{x), to distinguish it from the 
original free field, (po{x). Whereas ({>q{x) embodies the Uncertainty Principle, 
the infrared truncated free field does not, 

$o(a;),<l'o(x')] =0. (41) 

Note also that, whereas the VEV of v^o(^) diverges, the VEV of $o(^) "^^es 
not. We have proved that these very different quantities nevertheless produce 
precisely the same infrared logarithms. 

We have already seen that counterterms can be dropped at leading log 
order. The resulting, simplified Yang-Feldman equation is accordingly, 

$(t, x) = $o(t, ^) - ^ ll dt' V'($(t', x)) . (42) 

Note that we have called the field $(a;) to distinguish it from ^{x). This 
field $(x) commutes with $(x') for all x'^, just like the free field, $0(2^)- 
Further, VEV's involving it are competely free of ultraviolet divergences. 
They nevertheless agree exactly with VEV's of ^p{x) at leading log order. 

Taking the time derivative of (l42l) gives Starobinskii's Langevin equation 
[39], 

$(t,f) = $o(t,^) - ^y'{^(t,x)) . (43) 

Starobinskii's stochastic noise term is the time derivative of the infrared 
truncated free field, 

$o(t,x) = j J^^Ha - k)-={e^'Mk) + e-^^-at(fc)} . (44) 



(27r)3 ^ 

A simple calculation reveals that it behaves like white noise, 

n 



<l>o(t, f )$o(t', x) n) = —6{t - f) . (45) 
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2.5 Nonperturbative Solution 



Langevin equations of the form (1431) have been much studied [19]. Expecta- 
tion values of functionals of the stochastic field can be computed in terms of 
a probability density g{t, (p) as follows, 



Fmt,x)] 



+ 00 



d4>g{t,(P)F{4>) . 



(46) 



The probability density satisfies a Fokker-Planck equation whose first term 
is given by the interaction in fHSl) and whose second term is fixed by the 
normalization of the white noise (HSl): 



3H 



1 ^2 



^2 06^ 



47r2 



g{t, 



(47) 



To recover the nonperturbative late time solution of Starobinskii and Yoko- 
yama ^39] one makes the ansatz, 



lim g{t,(j)) = poo(0) 



(48) 



because the —V'{ip) force should eventually balance the tendency of infla- 
tionary particle production to push the scalar up its potential. This ansatz 
results in a first order equation, 



dgoo{(f)) 

Qoo{4>) 

The solution is straightforward, 

Qoo{4>) = Nexp 



2,m 



V'{(j))d(j) . 



(49) 



v{4>) 



(50) 



3 Generalizing Starobinskii's Formalism 

In this section we discuss the problem of generalizing Starobinskii's formalism 
beyond scalar potential models. We begin by explaining the two complicat- 
ing features: derivative interactions and passive fields which do not produce 
infrared logarithms. We then derive a formula which gives the leading log- 
arithm series for a general interaction. By applying this formula to simple 
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models it emerges that ultraviolet divergences can contaminate even the lead- 
ing infrared logarithms. This is because passive fields — and difi^erentiated 
fields of any type — make contributions of order one which multiply the lead- 
ing logarithms contributed by other fields. Unlike the infrared logarithms, 
these order one contributions derive from all portions of the free field mode 
sum, and from the full free field mode functions. Hence it is not valid to 
stochastically simplify passive fields as we did in the previous section. The 
correct procedure is to integrate them out, and then stochastically simplify 
the resulting effective field equations. We show that this amounts to com- 
puting the effective potential. 



3.1 Active, Passive and Differentiated Fields 

We use the term active to denote a field whose mode functions are right to 
produce infrared logarithms. Fields whose mode functions cannot produce 
infrared logarithms are called passive. A typical passive field is the massless, 
conformally coupled scalar. 



(51) 



The plane wave mode functions for this field can be worked out for any scale 
factor. 



v{t, k) 



exp 



-ik 



* dt' 



a{t') 



(52) 



The shall henceforth make the speciahzation to de Sitter so that. 



de Sitter 



v{t,k) 



1-— 
a 2 



exp 



ik ik -\ 



(53) 



These mode functions are not singular enough for the expectation value 
of two free fields to produce an infrared logarithm for D > 2, 



d'^-'k 
(27r)^-i 

(ggQ^-T 
(47r) 2 



„ife-A£ 



v{t,k)v*{t',k) 



(54) 



ik ik 







°° ^ ^ (kAx)eHa flo^ 
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Nor, again for D > 2, can one get an infrared logarithm from the vertex 
integration of the retarded conformal Green's function, 



cf (-^J 



H Jo 



dt' 



^ a ^ ^ a ^ 



(56) 



The absence of infrared logarithms from either source for the conformally 
coupled scalar is obviously related because its mode functions obey the same 
Wronskian as do the minimally coupled scalar mode functions. 



v{t,k)v*{t,k)-v{t,k)v*{t, k) 



-.D-l 



u{t, k)u*{t, k)-u{t, k)u*{t, k) . (57) 



This relation requires that a total of D — 1 factors of a must be shared 
between any two linearly independent solutions. For example, the real and 
imaginary parts obey. 



- 2Re(v)Im(v) + 2Re{v)lm{v) 



1 



-.D-l 



-2Re{u)lm{u) + 2Re{u)lm{u) . 

(58) 

For the VEV of a pair of free fields to produce an infrared logarithm requires 
that the far infrared {k ~ 0) mode function should approach a phase divided 
by k^~ , with no dependence upon a. If we make the phase zero then this 
fixes the small k dependence of the real part and we can use the Wronskian 
to infer how the imaginary part depends upon k and a, 



Re{u) 



lm{u) 



1 



k 2 



2#(D - 1)/^ a^-i 



(59) 



In this same limit the measure factor times retarded Green's function is the 
Fourier transform of the combination. 



u{t, k)u*{t',k) -u*{t, k)u{t',k) 



{D - l)H 



/a'\-D-i 



• (60) 



In contrast, the leading small k behavior of the conformally coupled mode 
functions is less singular than u{t,k) by a factor of {k/ Ha)~~^ . This pre- 
cludes getting an infrared logarithm from the VEV of a pair of free fields. It 
also shifts scale factors from the real part of v to the imaginary part. 



Re{v) 



; - - 



lm.{v) 



1 A;2 



(61) 
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Hence the measure factor times retarded Green's function goes to the Fourier 
transform of, 



v{t,k)v*{t',k) -v*{t,k)v{t',k) 



1 

H 



(62) 



As we saw in the previous section, positive powers of a' /a weight temporal 
vertex integrations overwhelmly at their upper hmits, which prevents them 
from producing infrared logarithms. Similar considerations apply to all pas- 
sive fields. 

Derivatives also suppress infrared logarithms. For if a free active field 
is differentiated with respect to space then its mode sum contains an extra 
factor of k and the logarithmic singularity evident in fl22|) is absent. A time 
derivative is even worse. In D = 4 it gives rise to two extra factors of k, 



u(t, k) 




ik 



ik 



H 



ik 



(63) 



Of course it is possible that a differentiated field in the Lagrangian will give 
rise to a differentiated vertex integration, rather than a differentiated free 
field, in the free field expansion. This also prevents the appearance of an 
extra infrared logarithm because the undifferentiated vertex integration is a 
function only of time at leading logarithm order. Hence a spatial derivative 
of it gives zero, whereas a time derivative would lower the number of infrared 
logarithms by one. 

Even though passive fields do not cause infrared logarithms, they can 
still transmit an infrared logarithm acquired through interaction with an 
undifferentiated active field. Fig. [1] depicts a two loop contribution to the 
VEV of Ff^iy{x)Fpf^{x) in which an infrared logarithm from the coincident 
scalar loop at x' is propagated through the virtual photon loop from x' to x. 

Passive fields can also induce interactions between active fields. For ex- 
ample, the photon loop in Fig. [2] induces an effective {ip*ipY interaction. The 
same comments apply to differentiated active fields. For example. Fig. [3] 
shows the differentiated scalar and the photon of the 3-point vertex fl78l) in- 
ducing an interaction between undifferentiated active fields. This is part of 
the full IPI 2-point function which has recently been computed at one loop 
order 
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X 




Figure 1: Two loop contribution to {Q\Fi^i^{x)Fp„{x)\Q) . 




Figure 2: Effective coupling in SQED. 




Lp*dLp dip*ip 
Figure 3: Effective ip*ip coupling in SQED. 
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3.2 Reaching Leading Logarithm Order 

One can understand the relation between coupling constants and leading 
infrared logarithms directly from the Lagrangian. Consider a term in the 
potential involving N undifferentiated, active scalars, V{(p) ~ cn^^ ■ An 
elementary exercise in diagram topology reveals that it requires two such 
vertices to add N — 2 loops to any diagram. Hence the loop counting pa- 

2 

rameter is {Cn)'^~^- Now consider the 2N extra fields associated with two 
A^-point vertices. In the free field expansion some of these 2N fields would 
contribute to retarded Green's functions while others would remain as free 
fields. It requires two free fields to produce a retarded Green's function whose 
associated vertex integration can result in an infrared logarithm. An infrared 
logarithm can also come from the VEV of a pair of free fields. It follows that 
the 2N extra fields from two vertices can produce at most A^ infrared loga- 
rithms. Hence the leading logarithm contributions to any VEV represent an 
expansion in powers of the parameter, 



Cn^" ^ [C'ijx\n'\a)\''-' . (64) 

For A(^^ the series would be in powers of A^ln^(a); for A(^^ we have already 
seen that the series is in powers of Aln^(a); for X^p^ the series would be in 
powers of A^ In^(a), and so on. 

Now consider a model which consists of an active field (p{x) and a pas- 
sive field tplx) that interact through a potential of the form, Kilj'''{d(p)"^ip"'. 
Because the interaction contains A^ = i+m+n fields, the addition of two ver- 
tices to any diagram increases the number of loops byA^ — 2 = i + m + n — 2. 
However, only the 2n active fields from these two vertices can contribute 
to infrared logarithms. Hence leading logarithm contributions to any VEV 
represent an expansion in the parameter, 

K^^difrip'' =^ [k'' X In" (a)] . (65) 

For example, the nhdhdh vertex of quantum gravity has i = 0, m = 2 and 
n = 1, which produces a series in powers of ln(a). The same result follows 
for any of the K^'h'^dhdh interactions of quantum gravity. 

A model whose leading logarithm solution has already been obtained is 
Yukawa theory [16], 
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+ (1 + 5Z2)ii^e^,i\d^ + '-A,,dr'')'ip^ - (/ + 5f)^ij^^g . (66) 

The scalar ip is active whereas the massless fermion is passive. The ba- 
sic interaction vertex is —f(pipipy/—g and comparison with (165!) shows that 
leading logarithm contributions represent an expansion in powers of /^ln(a). 
Contrast this with a —f^^^/—g vertex which would produce an expansion 
m powers of/2ln^(a). So we see that trading in an active field for a passive 
field, or for a differentiated active field, reduces the number of infrared loga- 
rithms per coupling constant. However, there will still be infrared logarithms 
as long as the vertex contains any undifferentiated active fields. 

3.3 The Role of the Ultraviolet 

It is instructive to examine the behavior of counterterms in Yukawa the- 
ory. Recall that they completely drop out at leading logarithm order when 
only undifferentiated active fields are present. The various counterterms of 
Yukawa theory behave as follows, 

6Zr^f , 6Z2^f , 5fr^f^ , , 5\^f\ (67) 

Field strength renormalization is irrelevant because the incorporation of ei- 
ther a. 5Z oi a. 5Z2 vertex would add a factor of without any undifferen- 
tiated active fields. The 3-point counterterm 5f is also subleading because 
it carries three extra factors of / with only one undifferentiated active field. 
However, the one loop contribution to the conformal counterterm 5^ adds a 
factor of with two undifferentiated active fields. These two fields could 
produce the extra ln(a) needed to remain at leading logarithm order. Simi- 
larly, the one loop contribution to the 4-point counterterm 5\ can also give 
leading logarithm corrections. Note that in neither case do we need to worry 
about higher loop counterterms because these would bring more factors of 

with no additional active fields. So we see that leading logarithm order in 
Yukawa theory requires one loop conformal and 4-point counterterms, but no 
other counterterms, and no renormalization at all beyond one loop. 

Passive fields engender ultraviolet divergences, even at leading logarithm 
order, precisely because they do not produce infrared logarithms. What they 
give instead is factors of order one which multiply the infrared logarithms 
from active fields. It is important to understand that these factors of order 
one derive from the full range of the free field mode sum and from the full 
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structure of the mode functions. Therefore, it is not possible to stochastically 
simplify Heisenberg operator equations which contain passive fields. The 
same considerations apply to differentiated active fields. 



3.4 Purging the Passive Fields 

It might seem as if there is no tractable formalism that describes the leading 
logarithm limit of models which contain passive fields. This is false. By 
integrating out the passive fields one obtains effective field equations which 
involve only active fields. Although these equations are hideously nonlocal, 
we will explain how they are equivalent to simple, local equations at leading 
logarithm order. 

Effective field equations are notoriously difficult on account of the field- 
dependent, inverse differential operators they involve. For example, consider 
integrating out a conformal scalar iplx) whose Lagrangian is. 



^ = -Id.^pd.ijg^-'V^ - - ^^V'v^ . (68) 

The resulting contribution to the ip{x) equation of motion takes the form. 



— Xip{x)(^x 



-g 



(69) 



where the quantity at the right is the coincidence limit of the ip propagator in 
an arbitrary ^{x) background0 We shall probably never know this Green's 
function for arbitrary ^{x). It can be expanded in terms of the conformal 
propagator iAcf (x; x') for = 0, 



X 



d^.iV^g^'d,) - 



1 / D-2 ' 



X 



/ d^x' 



iX\2 



iAc{(x;x) 



g{x')y^\x')[zA,,{x;x')] + (--) d^'x' ^ -g{x')v\x') 



X f d^x" ^-g{x") ^2(x")iAcf (x; x')iAcf (x ; x")i/\^i{x"; x) + . . . (70) 



"^It should properly be the coincidence limit of the ++ propagator of the Schwinger- 
Keldysh formalism [51] so that it depends only upon fields '^{x') in the past light-cone of 
x''. Giving this a proper explication would require a substantial digression which we shall 
forgo. 
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We do know iAc{{x; x') — in conformal coordinates it is just {aa'Y~^ times 
the massless, flat space propagator — but there is no way of doing the inte- 
grals and summing the series for arbitrary ^{x'). 

That we can obtain a tractable formahsm at leading logarithm order 
derives from two facts: 

1. To reach leading logarithm order, every active field must contribute to 
an infrared logarithm; and 

2. Integrating over passive field Green's functions cannot produce infrared 
logarithms. 

The first fact means we can ignore the spatial dependence of all the (/j's in 
expansion fl70|) . Of course their temporal dependence matters because this 
is the ultimate source of infrared logarithms. However, integrating infrared 
logarithms against a passive field Green's function does not change the lead- 
ing infrared logarithm. For example, compare the result of having infrared 
logarithms inside such an integral with the result of placing them outside. 



dt a 



d^'^x'GAx\x')W{a!^ 



HJo 



dt' 



In^(a) 



dt a 



I iD-l 



a 



In^(a) 
H 



dt' 



ln^(a' 

x'Gci{x]x') 



41n^fa 



D{D-2)H^ 



1 + 



ln(a) 



2 



41n^(a) 
D{D-2)H^ 



,(71) 



(72) 



To leading logarithm order there is no difference! 

lYiom. the preceding discussion we see that the leading logarithms are 
not changed by moving all the (/j's of expansion (ITOll from inside the various 
integrals to outside, 



X 



iX 



2 / iX 



X 



iArAx:x) 



--^^\x) / d^x'^)[iA,,{x-x')\ + (-y^^(x)) / d^x'^-g{x') 



X I d^x" ^ -g{x")i A^f{x; x')iAci{x'; x")iA^f{x"; x) + . . . . (73) 



21 



This last expression can be recognized as the coincidence hmit of the prop- 
agator of a conformally coupled scalar with mass rn? = ^ip^{x). A compact 

result for this can be given in terms of the parameter u = ^^Jl — 2\(p'^{x) / H"^ 



X 



I 



X 



H 



D-2 



■ D-l 



D-l 



D_ 

2 



(47r) 
H^-^ Til 



r(f) 



D 



f)r(i^ + .)r(i^ 



{An 



V) 



- - ^; i) , (74) 



(75) 



It remains to substitute fl75|) in (|69|) and recognize the result as minus 
the derivative of the unrenormalized, one loop effective potential, 



- KffM = -A<^(x) X 



H 



D-2 



Til - 



'D-l 



D-l 



(4vr: 



_D 

2 



(76) 



The factor of r(l — y) represents a one loop, ultraviolet divergence which 
can be absorbed by renormalizing the quartic conformal and quartic self- 
couplings. The model that results is of the scalar potential type ([8]) already 
solved by Starobinskii. 

A little reflection reveals that the technique we have just sketched must 
always work. Of course we can integrate out any subset of fields. The nonlo- 
cality engendered by doing this to passive fields will necessarily be restricted 
to passive field Green's functions. By definition of the field being passive, 
such Green's functions cannot produce infrared logarithms. Hence it is al- 
ways valid to extract undifferentiated active fields from inside the nonlocal 
effective field equations. The result must always give minus the derivative of 
the unrenormalized effective potential. 

A final point should be noted concerning expectation values of opera- 
tors, in the original theory, which involve passive fields. In computing the 
leading logarithm result for such an expectation value one must functionally 
integrate out the passive fields. This will produce a potentially divergent 
expression involving only active fields. The stochastic expectation value of 
this expression can then be computed using the effective potential. 
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4 Effective Potential of SQED 



The point of this paper is to obtain a leading logarithm solution for SQED 
and that is largely accomplished in this section. We begin by applying the 
procedure of section 3.2 to show that SQED gives a series in powers of ln(a) 
at leading logarithm order. We also work out which counterterms make lead- 
ing log contributions. We then integrate out the photon field and renormalize 
the scalar effective potential. The section closes with large field and small 
field expansions of the effective potential. 

We repeat the SQED Lagrangian iQ from the Introduction, 

£ = -il+6Z2)id^-teA^)ip*{d, + teA,)ipg^'''y^- 5^ip*ipR^ 
-^-{l + 6Z,)F,,F,,g^Pg''^^ - ^ (^^f ^ , (77) 

The complex scalar (f{x) is active whereas the photon is passive. The 

primitive 3-point interaction involves a passive field, a differentiated active 
field and an undifferentiated active field, so we have the case of ( IMl) with 

i = m = n = 1, 



ip*d^ifi - d^ip*ifi g'^'V^ =^ ln(a) . (78) 



The primitive 4-point interaction represents i = n = 2 with m = 0, so the 
result is the same, 

-e'^A^A^ip*^g'"'y^ =^ eHn{a) . (79) 

The various counterterms have the following dependences upon e^, 

6Z2 ~ , SZ3 ~ , 5^ ~ , (5A ~ (80) 

Hence field strength renormalization cannot contribute at leading logarithm 
order because the 6Z2 and 6Z3 interactions add a factor of with no undif- 
ferentiated active fields. On the other hand, the conformal and quartic coun- 
terterms correspond to £ = m = 0, with n = 2 and n = 4, respectively. We 
therefore conclude that leading logarithm SQED requires one loop conformal 
and 4-point counterterms, but no other counterterms and no renormalization 
at all beyond one loop. 

It is now time to integrate out the vector potential. Ever since the classic 
work of Coleman and Weinberg [SI] , it has been realized that this is is greatly 
facilitated in Lorentz gauge, 

dJV^g^''A^=0. (81) 
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The simplifications turn out to be even greater in de Sitter background so 
we will follow the usual practice, although other gauges can of course be 
employed ^55]. Dropping the field strength renormalizations, partially inte- 
grating, using the gauge condition and expanding the Lagrangian in powers 
of the vector potential gives. 



£ 

£o 
£i 
£2 



£0 + £1 + £2 



-d^(p*d„(pg 



SX 



2 ^ 



A,. 



-9 



(82) 
(83) 
(84) 
(85) 



Here is the vector d'Alembertian defined by \2^"Ai, = A^'P^. The effec- 
tive action is determined by the equation. 



[dA,]S d,{^g^^A, 



(86) 



Because £ is quadratic in the vector potential we can obtain the following 
explicit expression for F 

r[^*, ^] = 5ob*, Injdet v^(n'^^ - R^'' - 2eVV^7'^^) | 



d^xJ^ix) I d^yix 



fJ.U\ 



y)r{y). (87) 



Expression (1871) consists of three terms: the purely scalar parts of the bare 
action, a determinant factor, and the "current-current" term from completing 
the square in the exponential of the functional integral. One of the great 
things about Lorentz gauge is that this third term drops out at leading 
logarithm order. Recall from the previous section that leading logarithm 
contributions are unchanged by moving undifferentiated scalars around to 
different sides of inverse differential operators. But then a partial integration 
gives a gradient of Lorentz gauge propagator, which vanishes. For example, 
consider the leading logarithm contribution from the second term in J'^{x) / ie, 



ld''xip*{x)dpip{x)^-g{x)gP^{x){x 



-g[^^''' - R^"" - 2e V Vfi' 



y 



lead 
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—(p*Jd^xipdi 



-ggp^{x 



i 







y), (88) 



(89) 



-g [ni"' - R^"' - 2e V Vfl'^''] 

The same cancellation also occurs for the flat space effective potential [53]. 

The determinant term in (I57|) makes a leading logarithm contribution to 
the effective field equations, but only in the form of the coincidence limit of a 
massive photon propagator, with = 2e^ip*ip treated as if it were constant. 
The massive photon propagator and its coincidence limit have recently been 
worked out in de Sitter background [56] and the result is. 
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(47r)f lr(f + 1) 
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The parameter v has the following definition. 



D - 3x2 



D-3x2 2eVV 



(90) 



(91) 



(92) 



We can recognize the derivative of the unrenormalized effective potential 
from expression (11021) . 



K'fr(¥^V) = 6iD{D-l)H^ + yy^V 

fr(D-i) 



-D{D-l 



H 



D-2 



(4vr) 



r(f- 



-i: 



r 



-)■ 



-V) 



T{\ + u)V{\-u) 



(93) 



Renormalization is accomplished by first setting D 
the final term, 



rn? 



r(3- 



r(3-f) 



r(-2+; 



r(|+z.-f)r(|-z.-f) 



4 — e and expanding 



1 + 



r(| + ^)r(| 

eVV\r,/3 N /3 X 3 



+ 0(e) . (94) 
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The digamma function ipi^) has the following definition and expansions for 
small z and large z jlS] , 

^ ln(r(z)), 



dz 



-7 + E(-ircH(^-i; 



n-l 



n=2 



ln(.) - ^ 



+ 



12^2 120^4 

We make the following choices for the two relevant counterterms, 



5^ 



6X 



(An) 
{D- 



4 



^ + i + 0(D-4) 



(4vr)i 



4 



h 7 

D ^ 2 



0{D - 4) 



(95) 
(96) 

(97) 

(98) 
(99) 



Substituting (JHl) and (I98H99D in dHS]) and taking the limit e 
-2 3/7^ 



gives, 



87r2 



-1 + 27+ (-3 + 27 



+ 14 



'4- 



-V 



>.(100) 



Of course the divergent parts of 5^ and are fixed. Our choices for the 
finite parts are motivated to make the (v?*^:')^ and {}{>*^Y terms in ^^^(^999*) 
vanish, which keeps the scalar light as long as possible. Interestingly, the 
same choice for 5.^ cancels the leading infrared logarithm in the two loop 
expectation value of i{>* {x)i-{){x) p!3], and also keeps the one loop scalar mode 
functions from receiving any significant late time correction |57j . 

It remains to work out the effective potential and expand it for large and 
small field strengths. From fllUUp and the definition fl^2]) of v we see that the 
result depends upon the combination. 



~1p~ 



(101) 



Integrating (llOOp gives 
3//^ { , 



eff 



(_l + 27);2+(--+7)^2 

+ jdx{l+x) + |. (102) 
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An explicit power series expansion can be obtained for Ves in terms of the 
parameter, 



Az 



VI -8z 

2 2 



Substituting 



m 



2z + Az"^ + 16z^ + 0{z^) 
and performing the integral gives, 



(103) 



eflf 



87r2 
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2 ^ ^24 12 8 
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= \-- -C(3) 



^ 2m + 2 
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|A^2m+3 



2m + 3 2m + 4 



1-^C(3)]a^^ + 0(A^^) 



(104) 
(105) 



As already stated, our choices for the finite parts of 6^ and 6X cancel the 
order z and z"^ terms in the small field expansion. 

The large field expansion derives from substituting fl97j) in fll02p . 



3H^ 



eflf 



1 r 7 1 1 

-z^ ln(2; + 1) + [-^ + 2 ^^(2) + 7J + ^ ln(2; + 1) 



+ 



13 n 19 1 

-— + ln(2) + 27] ^ + — ln(z + 1) + 0(1) . (106) 



Because z grows for small H, as well as for large (p*(p, the {ip*(p)'^ln[ip*ip) 
term in fll06p should also agree with the flat space result of Coleman and 
Weinberg [SH \5E\. From equation (4.5) of their paper we see that it does. 
The full asymptotic expansion of (11061) up to order one also gives a generally 
accurate approximation of the potential, even for 2; < 1, as can be seen from 
Fig. H and Fig. [51 



5 SQED Stress Tensor 

The stress tensor of SQED is a composite operator which involves the pas- 
sive field Afj^. One must therefore integrate out and simplify the resulting 
functional of (p* and ip before its expectation value can be computed stochas- 
tically. That is the task of this section. As a bonus we obtain independent 
results for the two gauge invariant operators which principally comprise T^j^: 
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Exact potential (lines) and its asymptotic form (crosses) 




00 



Z 



the field strength bihnear and the scalar kinetic bilinear. The section closes 
with a discussion of the curious fact that the simplified stress tensor is not 

quite -g^^uVesi'^*'^)- 

Recall from the previous section that field strength renormalization coun- 
terterms do not contribute at leading logarithm order whereas the conformal 
and quartic counterterms do. The relevant part of the SQED stress tensor is 
therefore, 



T, 



ylv^V) V , (107) 

where the covariant derivative is D^^Lp = (d^ip + ieA^ip). In de Sitter back- 
ground the Ricci tensor is R^i, = {D — l)H'^g^y, and differentiated scalars 
with the same power of as undifferentiated scalars are guaranteed to be 
subleading logarithm. Hence we can simplify the counterterms to, 



'^K[^*^{R^lv-^g^,uR) +g,,u{v*vy% - iv*^) 



5\ 



^D-l){D-2)5iH^V*V+ —{V*^f\g,, . (108) 

It remains to integrate the vector potential out of the field strength and scalar 
kinetic bilinears, 

F^pFp, and {D.ipYD.ip. (109) 

Consider the general case of integrating the vector potential out of some 
operator 0[ip* , ip, A\ to obtain a new operator V^] depending only upon 

the scalar. 



[dAMd^{V^g^^A,)]e'^^^'^^^^WO[ip\^,A] = (^^^^'^^WO[ip\^] . (110) 



The functional integration is trivial because the Lagrangian of SQED is 
quadratic in the vector potential. After some partial integrations and ap- 
plications of the Lorentz gauge condition flHTl) . it consists of a part Cq which 
depends only upon the scalar, a linear part of the form Ci = J'^A^, and a 
quadratic part £2 = \A^V^^" Ay. We can read off the current from ([HI 
and the differential operator from ( l85i) . 
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ie df,if*if-!f*df,(f g 
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One evaluates the functional integral flllUp by completing the square, 







A^ + 



-r 



2 Vf^" 



;ii3) 



The operator 0[(p* ^ ip\ will therefore be the sum of terms from 0[ip* , if, A\ in 
which all combinations of the following replacements are made, 



A^{x) 
A^{x)A^{x') 



d^x' (x x'^J^{x') 



X 



X 



;ii4) 
;ii5) 



As we saw from expressions fl88tl89p of the previous section, the Lorentz gauge 
condition means that there is never a leading logarithm contribution from 
replacement flll4p . Hence we need only consider replacement flllSp . 

At this stage we must digress to discuss the photon propagator. We shall 
never know the inverse of the differential operator flll2p for arbitrary y:>{x). 
However, leading logarithm results only involve this operator evaluated for 
the special case where 2e'^(f*{x)!f{x) is a constant we shall call m^. That 
Green's function we do know [56j. It can be expressed in terms of the de 
Sitter invariant function of conformal coordinates = {f],x), 
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\x — X 
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The massive, Lorentz gauge photon propagator takes the form 



X 



B{y) 



dxt^dx'' 



+ C{y) 



dy dy 
^x^^ dx'" 
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The functions B[y) and Ciy) can be expressed in terms of a single function 

7(2/), 



B{y) = 
C{y) = 



1 



A{D-l)H'^ 
1 



A{D- 

The function 7(7/) is. 



-{Ay-y^)i{y) - {D-l){2-y)^{y) 
[2-y)i{y)-{D-l)^{y) 
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D + l D + l D ^ ^ 
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(120) 
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Here the parameter v is, 



V 



;i2i) 



Because our work is limited to coincidence limits we require only the 
integer powers in the Laurent expansion of 7(2/), 



-) 



m 



V{n + D-l)fy 



2 n=0 

2± 

2 



r(n + f + 



r(f -l)r( 2-f ) V[n + !^ + v)V{n + !^-v) ^y^n 
2 ' -/-va 



rfn- 



D 



..(122) 



When x'^ = x'^, the function y{x] x') vanishes. Because dimensional regular- 
ization ignores all D-dependent powers of zero, the coincidence limits of 7(2/) 
and its derivatives derive from factors of ?/° = 1. The two we require are. 



7(0) 



7'(0) 



D-i.H^ H^-'' (r{D~i) r(-f)r(^+z/)r 



•D+l 



-) 
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And the two coincidence limits we need of the photon propagator are. 
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(125) 
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D + l 
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7'(0)+7(0) 



D-1 



7'(0) 



D-1 



7'(0)-7(0) g^^g,,. (126) 



Here Z)p is the covariant derivative operator defined by DpA^ = A^-p. 

We can now integrate the vector potential out of the field strength and 
scalar kinetic bilinears. First, recall that the ordinary derivatives in the field 
strength tensor can be replaced by covariant derivatives. 



D^A, 



D,A^ 



(127) 



Combining this with relations (11261) and (1123111241) . and with our earlier in- 
sight about the J" terms dropping, we see that the leading logarithm form 
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of the field strength bihnear is, 



ein-fi*,v^] X lim \ DaD' 
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The analogous result for the scalar kinetic bilinear is, 



g-«r[(p*,(/p] ^ 
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r(i+z/)r(i-z/) 



(131) 
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Before computing the stress tensor we should comment on the explicit 
perturbative computations which have been done to check the field strength 
bilinear (11301) and the scalar kinetic bilinear (11331) . Of course there is no way 
to check the nonperturbative information these expressions contain! How- 
ever, we can compare against explicit one and two loop computations by 
expanding the parameter v in powers of m? = 2e^ip*(p, 



D-3 



1 2 4 

1 m ^/m \ 



D-3H'^ 



(134) 

From (11301) we see that the leading logarithm result for the field strength 
bilinear should be, 



Fc,i3{x)Fp„{x 
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T{D-1) 
(4vr)f r(f + l) 
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,(135) 
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+V^(2) + ^(^-l) - ^(2-:^) (Az.) + (0(Az.2)) [.(136) 



Now substitute the leading stochastic result using relation ([2 
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Because the two loop result is divergent we report only the divergent part, 



lead log 



167r2 



1 e 
X 

D-A TT 



^ln(a) + 0(e^ln'(a)) I. (139) 



The one loop (order one) result is trivial and is not, in any case, a check 
of the technique because it contains no infrared logarithm. The two loop 
(order e^) result agrees with the diagram given in Fig. [1] but we have not yet 
proved that the other two loop diagrams fail to contribute divergent leading 
logarithms. 

The scalar kinetic bilinear has been more thoroughly checked. Our sto- 
chastic prediction for it is. 



X • 



lead log 

r(i^-i) 

r(f+i) 



9pa 



D-lx H 



D 



4 



D 
2 



2 



(4vr; 



f r(#+r 



(47r) 

r(D-i-Az/)r(2+Az/) 
r(f-i-Az/)r(2-f +Az/) 

i,{D-l) 



.,(140) 



-m - ^(^-l) + ^(2-^)1 (Az.) + (0(Az.^)) 



;i4i) 



34 



2 TT2D-4: 
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1) + ^/'(Z^-l) - V(2) ln(a) + 0(eMn2(a)). (142) 



This agrees exactly, and for arbitrary dimension D, with the infrared loga- 
rithm in equation (146) of our recent, explicit two loop computation of the 
scalar kinetic bilinear [T3]. 

We can now assemble the various constituents of the stress tensor. From 
(11081) and we that the conformal counterterm contributes. 



(D-l)(Z}-2)5ei^W^7M- = X {D-\\D-2) 
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(143) 

Recall that z = e'^ip*ip/H'^ and e = A — D . We shall keep the same form 
as fll43p for the each of the four terms. From (11081) and fl^ the 
counterterm gives. 
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The field strength contribution results from combining f ll07p with (11301) . and 
it turns out to be finite, 
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The contribution of the scalar kinetic term comes from substituting (11331) in 
(11071) and then making use of the same expansion flM|) as for the effective 
potential. 
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Of course the divergences in (11431) and (11441) cancel those in (I150p , at which 
point we can take D = A. The final result has the form —g^iyVs{z) where, 
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It is apparent from Fig. that the stress tensor potential Vs{(p*(p) does 
not quite agree with the effective potential Ves{'^*'^)- This same sort of dis- 
agreement was also noted in the recent leading logarithm solution of Yukawa 
theory |16]. In both models the difference arises because the two poten- 
tials describe different physical processes: V^s controls the scalar's evolution 
whereas Vs controls the gravitational back-reaction. The two are distinct be- 
cause almost all the factors of in Ves are i?/12 for a general metric, and 
this changes the stress tensor, even in de Sitter background and at leading 
logarithm order. 

To see the point, consider a contribution to the matter Lagrangian of the 
form, 

AC = -F{R)y^ . (152) 
The corresponding contribution to the stress tensor is, 

' (163) 



AT, 



= 2R^,F'{R)-g^,F{R) + 2g^,F{RY%-2F{R),^,. (154) 

In D = 4 de Sitter background = g^j^^R/A, and we can ignore the deriva- 
tive terms at leading logarithm order. 



AT, 



F{R) - \rF\R) 



(155) 
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If all the factors of in V^g were R/12 for a general metric, then for de 
Sitter background we would have, 



F{R) 



87r2 V12 



R 



F{R)--RF'{R) 



= — -x-zf (z) . 

(156) 

The actual relation between Ves and V"s is tantalizingly close to ( 11561) . If 
we extract a factor of Si/^/Svr^ from each potential. 



3H^ 



and 



/s(^) , 



(157) 



then comparison of expressions (11021) and (I15ip reveals the following relation 
between the two dimensionless functions of z. 



1 



1 



fs{z) = 7^zf\z) - -z- -z 



(158) 



The reason for the extra contribution of —z/2 — 2:^/4 is that a small portion 
of the dependence in Kg is really the constant A/3, rather than R/12. 
These are the finite factors of In(if^) which derive from our counterterms 
through the expansion. 
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A-D 



-^-\hh') + o{a-d) 



(159) 



If we regard the arbitrary generalization of these factors of In(if^) as ln(A/3), 
rather than ln(i?/12), it corresponds to adding the following term to F{R), 



AF{R) 



3 /i?N2 



12eVV l/12eVV\2 
^ 2V 
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Because the logarithm vanishes for R = 12H^, the only change in F — RF' /2 
in de Sitter background is precisely the required deficit term. 



AF(i?) - -RAF'{R) 



dS 



?>H^ r I 1 21 
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6 Nonperturbative Predictions 

It would be silly to stop without exploiting the formalism we have developed 
to answer nonperturbative questions about SQED. One would like to know: 
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1. How large does the scalar field strength become? 

2. What is the asymptotic late time value of the photon mass? 

3. Does the scalar remain light? 

4. Does the vacuum energy increase or decrease, and by how much? 

5. What becomes of the electric and magnetic field strengths? 

Answering these questions is the task of this section. We begin by making 
the trivial generalization of Starobinskii's formalism from a real scalar to a 
complex one. We then exploit the results of the previous two sections to 
compute explicit answers to each of the five questions. 

The stochastic formalism of subsections 2.4 and 2.5 has a straightforward 
generalization to a complex scalar. One simply decomposes the complex field 
into two real scalars in the usual way, 

^{x) = + ^M^)) ■ (162) 

Now suppose expectation values of the quantum fields ipi{x) agree, at leading 
logarithm order, with those of the stochastic random variables ^i{x) which 
obey the Langevin equations, 

^^ = ^^0-—^. (163) 

The fields $jo(^,^) are independent sources of Gaussian white noise, 

$,o(i, x)^j^{t', x)) = ^5{t - t%j . (164) 

Then the expectation value of any function of the <l>j is given in terms of a 
probability density g{t, 01, ^2), 

d(f)J d(f)2git,(f)i,(f)2)F{(f)i,(f)2) . (165) 

-00 J —00 

This probability density obeys the Fokker-Planck equation, 

. 1 >f d r^Kff 1 1^ fH^ X ..... 
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We now exploit the fact that the effective potential depends upon the ipi 
only through the combination, 



li^l + ^l)- (167) 



2 

It follows that the probability density has the same form and that we can 
write the Fokker-Planck equation as, 

Here a prime denotes differentiation with respect to the variable 0*0 = (0f + 
0i)/2. 

It is apparent from the large field expansion (11061) . and from Fig. IH that 
Ves is bounded below. We can therefore make the ansatz of Starobinskii and 
Yokoyama [39] that the probability density approaches a time independent 
form at late times, 

lim ^(t, 01,02) = ^oo (0*0) • (169) 

t—^oo 

Substituting this in our Fokeer-Planck equation fll68p and making a few 
simple inferences implies, 

gU<P*4>)VM*<P) = -^Q'ooi^*^) =^ f'oo(0*0) = iVe-|i^^-(<^*<^) . 

(170) 

The asymptotic probability density can be more simply expressed in terms 
of the function f{z) introduced in equation (11571) . 

^?oo(0*0) = iVe-^(^) , (171) 

where z = e^0*0/if^ and. 



2 

+ fdxil + x) 



(172) 



Hence the late time limit of any function of the operator ip*{x)(f{x) can be 
reduced to an ordinary integral. 



/OO /* oo 

d(f)^ j^d(t>2 F(0*0) f?oo(0*0) , (173) 

/"OO J-T'^ "y 

= 2'kN dzzF{—-)e-f^'K (174) 
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Operator 


Expectation Value 




1.6495 X /e 


VP 


3.3213 X n /e 




7.6308 X //e/e^ 




3.2991 X H"^ 


Ml ^ 


.8961 X 3e2/7787r2 




.7223 X 3H^/8tt^ 




-.6551 X 3H^/87i^ 


( -^/ii^ -^pcr ) fin 


-9.5246 X H^/Stt'^ {g^ipQua - QpaQvp) 



Table 1: Late time limits of expectation values of some important operators. 



It is of course impossible to obtain analytic expressions for integrals of the 
form (11741) with the function f{z) in (11721) . However, it is nothing these days 
to evaluate such integrals numerically. We have done this using the "Nlnte- 
grate" function of Mathematica [SH]. The normalization factor is about, 



2t:N 







^ (175) 



2.16603 

The scalar reaches a nonperturbatively large field strength, 

lim/^*(x)^(x)\ = ^x2nN dzz^e^^^'^ ^ ^ x . (176) 

i-oo\^ V V g2 7q g2 2.16603 ^ ^ 

The photon mass-squared is = 2e'^ip*ip, so (I176P means that it reaches 
the asymptotic value, 

lim ^ 3.2991 x . (177) 



This is explicit, nonperturbative confirmation of the conjecture by Davis, 
Dimopoulos, Prokopec and Tornkvist [l2l HS] that inflation induces a nonzero 
mass photon mass. Indeed, the nonperturbative result is about a hundred 
times larger than one loop computations [60], [HI |12|, EH [62l |63] . 

Because inflationary particle production would be quenched if the scalar 
were to develop a mass comparable to the Hubble parameter, it is impor- 
tant to check that the scalar remains light. The scalar mass-squared is the 
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derivative of the effective potential, 



Ml = Kff(v'V) • (178) 
Our result for its asymptotic value is, 

Therefore the scalar is always light compared to the Hubble scale, and the 
approximation of treating it stochastically with the massless mode functions 
is justified. 

Table [1] summarizes our results for late time limits of various operators. 
One operator of particular interest is the stress tensor potential, Vs{(p*(p), 
given in equation fllSip . In view of relations (1157111581) we can express its late 
time limit as, 

87r2 2.16603 ^ ' 

That the surprising sign is correct can be seen from Fig. [7] which gives an 
expanded view of the stress tensor potential. It should be noted that the 
sign is due to the two negative terms in the z integrand. With just the first 
term the result would be, 

X 27riV rdz - f (^)e-^(^) = +^ . (182) 



87r2 Jo 2 ' ' ' 87r2 

The physical interpretation may be that inflationary particle production po- 
larizes the vacuum, which lowers the energy of a charged particle in the 
medium provided the charge density is not too large. 

Note from Table [1] that the expectation value of the effective potential 
is positive. As explained at the end of the previous section, this is no con- 
tradiction with our result for because the two potentials answer slightly 
different physical questions. V^g controls the evolution of the scalar field 
strength whereas controls the gravitational response. In particular, our 
result for implies that inflationary particle production induces a small 
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stress tensor potential 




CO 



0.5 1 1.5 2 2.5 3 



fractional reduction of the asymptotic expansion rate 
3Hl ^ 3H^ 



SttG X — — X .6551 



87r2 



3H^[l - .6551 X . 



183) 
;i84) 



This is insignificant even for the highest scale inflation {GH < 10^ ) con- 
sistent with the normalized CMB quadrupole and with the current upper 
bound on the scalar-to-tensor ratio [64j. The sign of the effect is nevertheless 
intriguing. 

Because the field strength bilinear involves coincident passive fields it 
requires renormalization even at leading logarithm order. This is evident 
from the factor of r(2 — y) in (11301) . As might be expected, the divergence 
can be absorbed with terms proportional to and {ip*(p)'^, 
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With this definition of the divergent part, the remaining finite part has the 
following D = 4 limit. 
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\^ pv^ per 



fin 



-z - zil + z) 
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- -Vl-8z) + + I [gppg^., - g^^g.p) .(186) 



The late time limit of this quantity is strongly negative. 



lim((F,.F,. 
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20.6305 
2.16603 



[g^ipgucT g pa-gup 



(187) 



Of course the quadratic and quartic parts of fll86p can be adjusted with 
the renormalization condition but the negative sign of the large field limiting 
form is unambiguous. 



(f F 

\^ pv^ pa 



fin 



' -z{z + l)\Yi{2z + 2)--z- 



-0(1)] (, 



ua gpagup 



(188) 

This indicates that the inflationary production of charged scalars increases 
the electric field strength (for example, fi = p = and u = a = i), relative to 
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its vacuum value, while the magnetic field strength (for example, ^ = p = i 
and u = a = j) is decreased. This makes good physical sense. Although the 
average charge is zero, there is about about one infrared scalar in each Hubble 
value. The local electric field is necessarily dominated by this charge, so the 
square of the electric field strength should increase. On the other hand, the 
nonzero scalar field strength engenders a positive photon mass which drives 
down the magnetic field by the Meissner effect. From expression (I147p we 
see that the net electromagnetic contribution to the stress tensor is that of 
a negative cosmological constant at leading logarithm order. This indicates 
that the damping of vacuum fiuctuations in the magnetic field is greater than 
the enhancement of the electric field. 

7 Epilogue 

Infrared logarithms are the manifestation of enhanced quantum effects me- 
diated by massless, minimally coupled scalars and gravitons. We call these 
active fields. The continued growth of infrared logarithms must eventually 
overwhelm even the smallest loop-counting parameter. At this point pertur- 
bation theory breaks down and one must employ some sort of nonperturbative 
technique to follow the subsequent evolution. A reasonable approach is to 
sum the series of leading infrared logarithms. Starobinskii has developed a 
simple stochastic formalism which accomplishes this for any model of purely 
active fields with nonderivative interactions [SOj [39] . 

More general models possess two sorts of complications: derivative inter- 
actions and couplings with passive fields, that is, fields which cannot cause 
infrared logarithms. We still do not have a general technique for handling 
derivative interactions. One deals with passive fields by integrating them 
out and then stochastically simplifying the resulting effective action of active 
fields. This amounts to computing the effective potential. The theory then 
reduces to the form that Starobinskii has already solved. This reduction was 
previously accomplished for Yukawa theory [16j, and we have done it here 
for SQED. 

Note that one must integrate out passive fields from the VEV of any op- 
erator. This can result in ultraviolet divergences even at leading log order, 
as we found for both the field strength and the scalar kinetic bilinears. The 
reason for this is that passive fields contribute factors of order one which 
multiply the infrared logarithms contributed by active fields. Whereas in- 
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frared logarithms derive entirely from the long wavelength part of the free 
field mode sum, the factors of order one come as much from the ultraviolet 
as from the infrared. 

The nontrivial role of the ultraviolet also shows up in the fact that certain 
renormalization counterterms can make leading order contributions. In both 
Yukawa theory pB] and SQED the one loop conformal and quartic coun- 
terterms contribute. However, no other counterterms contribute at leading 
logarithm order, nor do any higher loop counterterms matter. 

We have obtained nonperturbative results for the VEV's of 
and Tfj_y. Table [1] summarizes these. Our result is that the scalar approaches 
a nonperturbatively large field strength. This confirms the conjecture of 
Davis, Dimopoulos, Prokopec and Tornkvist |^ The scalar remains 
perturbatively light, which means the computation is self-consistent. 

Our result for the stress tensor is curious in two ways. First, although 
it takes the form T^^ — > —QfiuVs, the stress tensor potential is not quite 
equal to the effective potential. This does not mean one potential is "right" 
and the other "wrong." Rather, they are both the correct answers to slightly 
different physical questions. The effective potential controls how the scalar 
evolves at leading logarithm order, whereas the stress tensor potential de- 
scribes how this evolution serves as a source for gravity. 

As shown in section 5, the difference between and V^fj arises because 
almost all the factors of in are actually R/12 for a general metric. 
One consequence is that V^g represents a peculiar modified gravity theory 
which may have important implications for cosmology. Lagrangians of the 
form F{R) seem to be theoretically viable, and they can be tuned to give any 
desired evolution for the scale factor [65j. However, there is no justification 
for such models from fundamental theory. In contrast, the modified gravity 
model we get is uniquely fixed and thoroughly justified — although it may 
not, of course, do anything interesting. 

The second peculiar thing about our result for the stress tensor is that 
it reduces the vacuum energy. The physical interpretation for this may be 
that the inflationary production of charged scalars polarizes the vacuum, 
which lowers the energy of charged particles in this medium provided the 
charge density is not too large. Supporting this conjecture is the fact that 
the electromagnetic contribution to the stress tensor is that of a negative 
cosmological constant. 
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